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Abstract: We make a first application of the linear gravitational potentials of the conformal gravity theory 

to the distance scale associated with clusters of galaxies, with the theory being found to give a reasonable 

value for the mean velocity of the virialized core of the typical Coma cluster, with no need to invoke the 

^ ' existence of dark matter. 

< 

In a recent series of papers (Mannheim 1994 gives a full bibliography) Mannheim and Kazanas have 
explored conformal invariant fourth order gravity as a candidate alternative to the standard second order 
Newton-Einstein theory. In their initial paper they showed that in this theory the potential of a gravitational 
' source takes the form V(r) = —(3c 2 /r + r yc 2 r/2, to thus both nicely recover the Newtonian potential and 

at the same time introduce a new source-dependent length scale 7 which would then parameterize any 
possible departures from the standard phenomenology. Moreover, Mannheim (1993) then showed that with 
the use of this potential it is possible to fit the rotation curves of some typical galaxies without needing to 
assume the existence of any dark matter, with the contribution of the linear potential when integrated over 
the observed surface brightness data then completely replacing the dark matter contribution in the same 
galaxies. Intriguingly it was found from the fitting that the strength of the linear potential of an entire 
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galaxy was typically of the order of the inverse Hubble length to suggest a possible cosmological origin for 
7. Armed with an actual magnitude for the galactic 7 we can now apply the conformal theory to the first 
available distance scale beyond galaxies, namely that of clusters of galaxies. Interestingly, we find that the 



conformal theory deviates from Newton there by just the amount needed to nicely accommodate a virialized 
cluster core without invoking dark matter, though the theory may turn out to have some difficulties should 
& [ entire clusters prove to be virialized. 

In applying the standard cluster virial analysis to theories with rising rather than falling potentials two 
problems arise which are not considered in the standard theory. First since the standard discussion generally 
considers the two-body collision dependent term in the Boltzmann equation to be a local perturbation on 
the global mean field set up by the gravitational field of the rest of the galaxies in the cluster we have to 
reexamine the entire formalism in light of the conformal gravity linear potentials which grow with distance 
and which can thus not be thought of as producing localizable collisions at all. Moreover, in a theory with 
rising potentials, we are not free to ignore the effects due to all of the rest of the galaxies in the Universe, 
so that even if a system such as a cluster of galaxies is geometrically isolated, that does not immediately 
mean that it is gravitationally isolated in our theory or that it is bound purely under its own self forces. 
As regards the inapplicability of the Boltzmann analysis within clusters in the linear potential case, we 
note that fortunately there is an altogether far more general statistical analysis beyond Boltzmann, namely 
that based on the Liouville equation as conveniently treated via the very general Bogoliubov, Born, Green, 
Kirkwood, and Yvon (BBGKY) heirarchy (a heirarchy which reduces to the Boltzmann equation only under 
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very special circumstances and which remains valid even when no such reduction is possible). Since the 
BBGKY hierarchy is valid for any central potential it is thus immediately valid in the linear potential case 
too, with our ability to apply the standard Jeans and Vlasov virial equations then reducing to the degree to 
which we can ignore two-particle BBGKY correlations (see Mannheim 1995 for details). Since the BBGKY 
relaxation time is currently unknown it is not immediately clear how much of a typical cluster has so far 
virialized (we argue below that possibly only the core is virialized), and so we develop below a formalism 
which can deal with both fully and partially virialized clusters. As regards the effect of the global coupling 
of the cluster to the rest of the Universe, this would require the developing of a theory for the growth of 
inhomogencitics and galaxy formation, with the relevant issue for motions within clusters then being not 
so much the coupling of each cluster to the general Hubble flow produced by a homogeneous background 
distribution of linear potential sources, but rather its coupling to the deviations from that flow caused by the 
presence of inhomogeneities. Since a theory for the growth of inhomogeneities in conformal gravity has yet to 
be developed, we are currently unable to address this issue explicitly, though as a first step we will determine 
approximately where it is that the cluster begins to merge with the general cosmological background, and 
should thus immediately anticipate that only those regions of the cluster which lie well within this merging 
radius are well enough isolated to be able (if then given enough time) to decouple and subsequently virialize. 

To treat clusters in general let us consider a spherically symmetric cluster with matter volume density 
a(r) and matter surface density I(R) contained within some volume of radius Rm, and let us suppose that 
only some inner kinematic region of the cluster within some virialization radius r m has so far had time to 
virialize. Then within the r < r m region we may ignore all two-body correlations with the Jeans equation 

|>(r) < vl > ) + ^l( <v 2 > - < vl >) = -a(r)V'(r) (1) 

then holding in this region. Despite the fact that Eq. (1) only involves points with r < r m , we note that the 
potential V(r) needed for it is obtained by integrating the Newtonian and linear potentials over the entire 
cluster and not just the virialized region, with the linear piece (unlike the Newton piece) actually receiving 
contributions from points with r > r m . Now observationally we measure the two-dimensional projected line 
of sight velocity distribution average < cr 2 (R) > which is related to three-dimensional velocity averages (the 
ones which appear in the Jeans equation) according to 

I(R) < 4(R) >= 2 £ M r(r f_ g g )1/2 (< vl > (r 2 - R 2 ) + < v 2 > R 2 ) (2) 

with the integration range all the way to the cluster radius Rm in Eq. (2) thus involving not only the desired 
virialized r < r m region but also the non-virialized r > r m region where Eq. (1) does not hold. However, 
integrating Eq. (2) itself over a sphere of radius r m then yields (Mannheim 1995) 

2ir / dRRI(R) < o*{R) >= — / drr 2 a(r) (< v 2 > +2 < v 2 g >) (3) 
Jo Jo 

an exact and completely general relation in which the undesired unvirialized r > r m region has crucially 
been projected out. Since Eq. (3) only involves the region r < r m , it thus only involves the region where we 
have presupposed the Jeans equation to be valid, with its use then yielding for the spatially averaged line of 
sight velocity 

Mv 2 p (r m ))av dRRI(R) = ^ drr 2 a{r)rV'{r) (4) 
to give finally the form of the virial theorem for partially virialized systems. 
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For the most studied cluster, the Coma cluster, the relevant data may be found in White et al (1993) 
and references therein. At the distance of Coma an arc minute is 20 jh kpc (for a Hubble parameter 
Hq = 100ft- km/sec/Mpc), so that the standard Abell radius is 75' for Coma. The surface brightness may be 
approximately fitted by a modified Hubble profile (~ 1/(R 2 + Rq)) with a core radius i? = 9.23' and number 
density normalization (Jq/Rq = 0.016 galaxies per cubic arc min., with the observed cluster data going out 
about to 3° ~ 20i?o = Rm or so where we shall thus cut off the distribution. (A modified Hubble has to be 
cut off somewhere since it would otherwise yield an infinite total mass). White et al (1993) quote a total 
blue surface luminosity within the Abell radius of Lb = 1-95 x 10 12 /h 2 Lbq, and a mean projected line of 
sight velocity of 970 km/sec for a convenient magnitude limited cut on the data which restricts to R < 120'. 
For such a mean velocity, the time required to cross the associated 240' diameter is 1.5 x 10 17 /h sec, which is 
of order 1/2H , i.e. of order half a Hubble time, and thus we should not expect the entire cluster to have yet 
had time to virialize. Hence, for the purposes of this study, we shall simply assume that only the inner region 
cluster core has so far virialized. Giving each galaxy an average blue luminosity of 5.99 x 10 9 /h 2 Lbq, then 
yields the requisite total 1.95 x I0 12 /h 2 Lbq surface blue luminosity within the Abell radius, to thus fully 
specify I(R). Using as typical the mass to light ratio M/Lb = 5.6IiMq/Lbq obtained for the galaxy NGC 
3198 in Mannheim(1993) enables us to determine the mass volume density associated with tr(r). It is very 
convenient to express this mass density in units of the standard critical density p c = 3Hq / 8ttG, and we find 
that cr(O') = 241.5p c , cr(56.8') = p c , er(120') = 0.11p c , and cr(185') = 0.03p c . The cluster is thus apparently 
merging with the general cosmological background at no more than 185' or so, and would be restricted to 
the first 57' if the density of the Universe is critical. Thus in a low density Universe we would put the edge 
of the cluster at 185', while in a high density one we would only consider the potentials of the first 57' of 
data as contributing to the velocity dispersion, with the next 128' of data then only contributing along with 
the rest of the galaxies in the Universe to the general Hubble flow. (Noting that the conventional estimation 
of the cosmological ratio p/ p c is made in comoving coordinates while our analysis here involves the same 
ratio in static coordinates, our determination of where the static cluster actually merges with the comoving 
background is thus perforce only a rough estimate.) Since the actual density of the Universe represents one 
of the key unknown issues in cosmology, we shall calculate cluster core virial velocities for both the high 
and low density Universe cases, and actually find below that the values that we then obtain turn out to be 
insensitive to where the cluster ends. Moreover, no matter where we put the edge of the cluster, we should 
not expect the matter within the cluster but close to this edge to necessarily be completely decoupled from 
the cosmological background, so that again only core virialization seems reasonable. However, since the 
modified Hubble is still a falling profile, it turns out that one quarter by volume of the entire 20i?o cluster is 
contained within r < 2.5i?o, while one half is contained within r < 5Rq. Thus virialization of only the core 
region is quite non-trivial. 

Taking V(r) = -f3 ga ic 2 /r + Jgaic 2 r/2 as the potential put out by a typical individual galaxy and taking 
TV to be the total number of galaxies contained in the entire cluster (and not just the number N(r m ) 
contained in the virialized r < r m region) fixes the overall normalization {Nf3 ga ic 2 /i?o) 1 ^ 2 of the Newtonian 
potential contribution to the cluster virial, while taking as typical the NGC 3198 gamma to light ratio of 
9.2 x 10^ 40 h 3 /cm/ Lbq also obtained in Mannhcim(1993) then enables us to fix the overall normalization 
(A r 7 g oiC 2 i?o) 1 ^ 2 of the linear potential contribution as well. Thus for a Coma cluster composed solely of 
luminous matter alone, the overall normalizations (iV7 ga ;c 2 i? ) 1 ^ 2 and {N(3 ga ic 2 /Rq) 1 ^ 2 needed for Eq. (4) 
take respective values of 10960 km/sec and 576 km/sec for a cluster cut off at Rm — 20R (N = 425 
galaxies). In the absence of any dark matter the luminous Newtonian contribution to the virial is thus 
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negligible, while the linear contribution associated with the luminous matter is substantial. Specifically, if 
the entire Rm = 2QR cluster is virialized, Eq. (4) yields a virial velocity a p (20R ) = 10178 km/sec, while 
also yielding partial virial velocities a p (Ro) — 1089 km/sec, a p (1.5Ro) — 1678 km/sec, a p (2Ro) = 2195 
km/sec, and <7 p (6.15i?o) = 5018 km/sec in various inner regions. Similarly, if we cut off the cluster at 
57' = 6.WR a (to yield N = 242 galaxies, (iV7 ffa; c 2 i? ) 1/2 = 8261 km/sec, (N (3 gal c 2 / i? ) 1/2 = 435 km/sec) 
we obtain the partial virial velocities <j p (Rq) = 1028 km/sec, a p (1.5Ro) = 1583 km/sec, a p (2R n ) = 2070 
km/sec, and <7 p (6.15i?o) = 4885 km/sec. The core region velocities are thus essentially insensitive to whether 
we use a high or low density Universe cut-off. From the line of sight velocity data points we find that the 
numerical average of all the points in R < 1.3i?o is 1200 ± 195 km/sec, while that of the all the points in 
R < 1.7i?o is 1185 ± 195 km/sec. However, before we assess the significance of these numbers, it is important 
to note that once less than the entire spherical cluster is virialized, then any given line of sight through the 
sphere, even those at small impact parameter R, will pass through both virialized and non-virializcd regions 
(the integration in Eq. (2) is to Rm and not merely to r m ), so that the detected projected velocity at that 
R will include some non-virialized contributions as well. For instance, if r < 2.5i?o is virialized, then out 
of a 20i?o cluster the percentage of line of sight material which involves unvirialized radii r > 2.5Rq is 25% 
at R — 1.5i?o, 44% at R = 2i? , and of course 100% at R = 2.5Rq- Thus even though the partial virial of 
Eq. (4) itself only involves integrating up to r m , the very use of initial input raw velocity data to estimate 
a magnitude for a virialized a p (r m ) in the r < r m region becomes suspect once the cluster is less than fully 
virialized. From our calculated virial velocities we see that conformal gravity would thus appear to have 
no difficulty accommodating a virialized inner cluster region of the order of r m ~ 1.5i?o without needing to 
invoke dark matter, and given the just noted limitation on the use of the input velocity data in partially 
virialized systems, the theory could possibly even accommodate up to r m <~ 2.5i?o, a region which contains 
close to one quarter by volume of all of the matter in the entire 185' of the cluster. Moreover, given the 
relevant time scales which were discussed above, it would even appear to be quite reasonable to expect inner 
region virialization up to one or two scale lengths or so. While we would certainly not expect any larger 
a portion of the cluster to have yet virialized, a first principles determination of the two-body correlation 
function and of its potential impact on Eqs. (1) and (4) could nonetheless prove to be very instructive, and 
might possibly even turn out to be definitive for the theory. (It is also possible to test the conformal theory 
in a way which is actually insensitive to how big a fraction of the cluster has in fact virialized, viz. cluster 
gravitational lensing which responds to all the matter in the cluster virialized or not; thus a yet to be made 
study of the conformal theory predictions for lensing should eventually provide an independent and definitive 
way of testing the theory on whole cluster scales.) Other than this issue though, it would appear that, in 
the first instance at least, the conformal gravity theory is indeed capable of meeting the demands of cluster 
virial velocity data, with the linear potential theory thus readily being extendable from galactic scales up to 
the much larger ones associated with the virialized core regions of clusters of galaxies without encountering 
any major difficulty. This work has been supported in part by the Department of Energy under grant No. 
DE-FG02-92ER40716.00. 
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